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ABSTRACT 


This work consists of the development of the six degree 
of freedom non linear model of a sea launched generic bank-to- 
turn cruise missile attacking a medium sized combatant ship. 
Two guidance and control schemes are compared in the terminal 
phase. The first, or baseline guidance scheme (pop out 
jmem@euver), uses a S0=fo0t altitude hold for an ingress phase, 
followed by a pop out maneuver, and then an attack phase which 
uses proportional navigation in elevation and azimuth planes 
along with bank-to-turn maneuvering. The second scheme 
(sea skimmer) uses identical ingress and attack phases but 
eliminates the pop out maneuver. Miss distances for both 
schemes are compared while varying missile roll rate limit, 


BOM blinking frequency, and burn through ranges. 
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IT. INTRODUCTION 


The relative merits of bank-to-turn versus skid-to-turn 
missiles have been argued for years. Bank-to-turn missiles 
which must roll to a commanded bank angle before a lateral 
acceleration can be commanded in the appropriate direction are 
inherently slower to respond to target maneuvering than 
skin-to-turn missiles of comparable lateral acceleration 
Bepoemiathes., This is particularly true if the bank-to-turn 
missile has limited roll authority or limited roli rate 
meeeeemey.s On the other hand, as discussed in Gonzalez, 
meter), tt 1S sometimes difficult to build a skid-to-turn 
Meesile which meets certain operational goals and performance 
Peeeemra and also has sufficient control authority in both 
the y and z directions to generate the large lateral accel- 
-eeerods required to perform skid-to-turn. 

In this thesis the terminal guidance accuracy of a 
bank-to-turn cruise missile is studied for a sea launched, 
Seamtarget cruise missile of conventional wing-tail config- 
Memonon. it 1S assumed that roll authority is extremely 
limited due to operational design restrictions which impose 
the use of low authority ailerons or differential stabilator 
(no ailerons) for rolling. Lateral acceleration is also 
assumed to be extremely limited in the y direction due to 
the presence of only a small vertical surface at the tail 
(vertical stabilizer and rudder). 


1% 





Two bank-to-turn guidance schemes are compared for 
accuracy against a medium sized combatant ship employing 
electronic countermeasures (an ECM blinker mounted aft on 
the ship). The first baseline scheme employs a pop-out 
maneuver consisting of a low altitude run in, azimuth offset, 
Pee-up and roil into the target using proportional navigation 
in azimuth and elevation. The second, or sea-skimmer scheme, 
Siploys a straight-in low altitude attack. Miss distances 
for both schemes are compared while varying missile roll 
Paoewtimio, CM blinking frequency and burn through ranges. 

oo eee Simulation was coded to perform the study. 
Missile dynamics were represented by 6 degree of freedom 
nonlinear equations of motion with table look-up aero- 
dynamic coefficients. These coefficients are representative 
of a generic cruise missile with a conventional wing- tail 
Seema cion Of limited roll authority. Inner loop 
augmentation and autopilot modes were designed to improve 
(mereebescamping in all axes, provide commanded load factor, 
commanded bank angle, and turn coordination. Autopiiot 
outer loops were designed to provide altitude hold, and 
vertical flight path angle hoid. Guidance loops were designed 
to provide proportional navigation in elevation and azimuth. 

pince this study involves the influence of flight 
dynamic parameters (roll rate limits) on terminal guidance 
accuracy, no extensive tracker modeling is employed in the 


Simulation. It is assumed that the seeker always tracks 
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baie point target perfectly. The point target is, however, 
moved by an ECM blinker and contaminated with a glint model. 

In Chapter II the equationsused to represent missile 
dynamics are presented along with the methodology used to 
Preece Up aerodynamic forces and moments. In Chapter III 
the development, design and simulation of the missile auto- 
pilot modes are presented. In Chapter IV the design of the 
guidance systems and a complete description of the two 
guidance schemes is presented. In Chapter V the CSMP 
Simulation is presented followed by the results and conclusions 


Seeoias suvudy. 
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Ti. DEVELOPMENT AND SIMULATION OF MISSILE DYNAMICS 


In this chapter the linear and nonlinear mathematical 
models are developed which are used to describe the flight 
dynamics of a generic bank-to-turn cruise missile. The 
linearized mathematical model of missile dynamics is used 
mmedesten the missile autopilot control laws and guidance 
laws. The nonlinear mathematical model is used to accurately 
represent missile dynamics in a CSMP computer simulation of 
a sea launched generic cruise missile attacking a medium 
eezea combatant Ship. Since the purpose of this investigation 
moembo perform a detailed evaluation of terminal control 
laws, it was decided that a linear dynamic simulation would 
MeemaaecWately represent missile motion; hence a full six 
degree of freedom (6 DOF) nonlinear dynamic simulation was 
pmeoded. the overall system that is modeled in this study 
meee eic pea in Olock diagram format in Figure 2-1 and will 


be developed in this thesis. 


A. MISSIL& BRQUATIONS OF MOTION 

The full nonlinear six degree of freedom rigid body, 
dynamic equations of motion are used to represent the motion 
Of a generic bank-to-turn cruise missile. The aerodynamic 
forces and moments are built up from representative tabular 


coefficient data for the generic missile. These tabular 


Z0 





data aregiven in Appendix A in graphical form along with the 


Po ueteme Micstle'’s phycsacal~and geometric characteristics. 
+t. Assumptions 

The assumptions used in the development of the 
equations of motion are given below. 

4. ‘The earth is flat, does not rotate, and is fixed 
in inertial space. 

TiomiictcsOr bone missile is constant. 

feitiew@imsstle 15 a4 rigid body. 

Gq. The mass distribution of the missile is constant. 

e. Engine angular momentum is neglected. 

f. The missile is symmetric about the body's xy 
plane. Therefore, the products of inertia Tey and tre, are 
zero. 

ait wclowmmMenbMrtey Line 1s parallel to the missile 
mire axis. thus, the thrust components rt and qT, a 
Zero. 

h. The density of the atmosphere is constant. 

ie ie wemeine thrust line passes through the missile 
Genter of mass. Therefore, the moments due to thrust Lins 
Min and Na are Zero. 

lemeeine wind is calm. 

k. Right aileron trailing edge down is positive 


deflection. Negative elevator deflection yields nose up 


pitching moment (positive). 
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2. Coordinate system 


An earth fixed coordinate system is established with 
meeormein fixed at the initial position of the missile. The 
Xe Paro wormos tO true north, the Le axis POHNLS east, and 
the Le axis points toward the center of the earth. Altitude 
(H) therefore equals the negative of Za. The system is 
considered to be an inertial system. 

A stability axis coordinate system is used with 
peemdidabes x, yY, 2 fixed at the missile center of mass and 
oriented such that the x axis lies along the vehicle's 
forward velocity vector in steady state balanced cruise 
Meesaiecusl,/> Mach number. The y axis is chosen perpendicular 
Wome plane Of Symmetry and 1s oriented out the right wing. 
iMtem2naxi1s 1S chosen perpendicular to the x axis in the down 
direction and in the plane of symmetry. 

Pe mom tguaboms of Motaon 

iow oMlovinoerwMondincar equations as developed by 

Hewett {Ref. 2} describe the motion of the missile in 


6 degrees of freedom. 


m (U-VR+WQ) = -mg sin 9 +X+T a.) 
m (V+UR-WP) = mg sin ¢ cos 9 + Y (1b) 
m (W-UQ+VP) = mg cos ® cos 6 + Z (1c) 
PI -(R+PQ)I, + RQ(To-yy) = by (2a) 
QI tPR(I,,-1,,) + (Pe RT) T = My (2) 
RI PIy,tPO(T,-T. ORI, = Ny (2c) 
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6 = Pt+t(Q sin 6 + R cos ©) tan 6 (3a) 
ao ouneoe= 9h sin 0 coos 
vY = (Q sin 6 + R cos 0) sec 4 (3c) 


Equations (1a), (1b) and (1c) describe the transla- 
Meme motion of the missile. Equations (2a), (2b), and 
(2c) represent the rotational motion. Equations (3a), 
(3b) and (3c) are the Euler relations for bank, pitch and 
yaw angles. 
Bouations (4a), (4b), and (4c) describe the missile's 


position referenced to the earth fixed system (Xo, Yas Bose 


EB 
. = U cos ¥eos 8+ Vicos¥siné@sind-sin¥ cos? ) (nen) 
+ Wlceos¥sinéOcosft+sinvsin® ) 
= = U sin Yeos 6+ V(sin¥siné@sind+cos¥cos® ) eb) 


+ W(sin¥sin@cos6-cosY¥sin® ) 


= Usin6-Vcos@siné-Wceoscosé (Ae) 


Je 


if 
Equations (5), (6) and (7) are also required. 
Equations (5) and (6) yield angle of attack and side slip 


wmoee respectively im terms of velocity components U, V and 


W. Equation (7) represents the total velocity of the missile. 


' 
« = arctan —— = arcsin a i (5) 
(U<+wW* )* 
eos wal a 9 = arcsin -~ Dee (6) 
U~+w (U +V +W ‘a 
Vn > (u*4+V~4W*) 2 (7) 


Figure 2-2 illustrates the positive directions and 
Piemuotemoteiorces, moments and velocities in the body 


coordinate system (x, y, 2). 


a4 





Aerodynamic and Thrust Moments Linear and Rotational (Angular 
Velocities 





Figure 2-2 Body Coordinate System Bie yeren pa. On 


£5 





Bquations (1a), (1b), (1c) and (2a), (2b), (2c) can be 


=~ 


rewritten in state format as follows: 


§ = -gsinO+VR-WQ + - + - (8a) 
Vo= gsin®cos9-UR+WP + - (8b) 
W = gcos®cos9+UQ-VP + - (8c) 
P + eee tT, -1,,)1 eg(i" <1... 4a / 
A 22 at XZ Dex ZL XZ ZZ yy 2% Xz 

ene (9a) 
= = My-PRUL, =a Iga) (PPR IT, i. (9b) 
R =| Boe Naty PQtT, I XX Peete) PRTG T Ty) | / 

(Th lag Tae) (9e) 


Now integrating equations (3), (4), (8) and (9) results 
in U, V, W, P, Q, R, % 8, ¥, Xa, Y, and Z, Sea wae 
Missile’s velocities, angular rates, and angles and positions 
have thus been described. 
4. Trim Rgquations 
Powe owcrromueane bevel flight at constant velocity 


the following variables are all zero: U, V, W, P, 9, R, 


om O >» Vy Lo» Pags Rags Rag» Pag a pee. tog Yoo? 
Bags Vags “ss Thus the trim equations reduce to: 
KogtT = mg sind,, 
Zea = -mMg cos Bac 
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Xm «= Vg g COBB gg COSB gg ~ Va gGSIN¥a gtW ag COS¥agGSING gg 


oc 
Tee UsgSintgs COSsdggtV gg cos¥agtW oosin¥osinGss. 
tanOgg = Wog/Ugg 


a Ps 
ee = KU age + OW De 
Tac $5 5S 


Since bank angle equals zero and the sum of flight 
path angle and angle of attack equal the pitch angle, then 
mmeverdy Stave pitch angle equals angle of attack. 


Ugg = “gg 


These equations are used to define the initial cruise condi- 
tion of the missile which is low level straight and level 
eruise flight at 0.75 Mach number. The precise initial 
mene CONdi1tLON and initial target position is described in 
Peeoemdix Db. 
peeeouner Userul Relations 
Components of load factor at the missile center of 


mass in the y and -z directions are given below: 


— ~Z a es a a 
oe - (Q } (10a) 
y Vo 
a E (8+ k) (100) 


ee caemeom of Oontrols and Control Limits 
ime control configuration considered in the thesis 
is a conventional wing and tail with conventional rudder 
(et), aileron (8) and elevator (n) control surfaces. In 


Picea Go) COver 4 missile configuration which uses 
Ba 





eer erentia! horizontal stabilizer (stabilator) for roll 
meron as Opposed tO aileron, the quantities left stabilator 
deflection (ny; ) MGmominrestapalator deflection (np) are 
defined and calculated. 

Limits are imposed on rudder, elevator, and aileron 
travel. Limits are not imposed on differential stabilator 
travel but stabilator travel is presented for all simulations 
to show the stabilator travel required if indeed that were 
the control configuration. The aileron power coefficient 
data presented in Appendix A reflects the limited control 
Peer iby in roll typical of a missile configured with 
meeeerenplal Stabilator for roll control as opposed to 
aileron. 

The following control definitions and limits are 
pee oe ed ; 

a) Thrust (T) is a constant and is oriented along 
five X axis. 

T = Ty smecens Lant Ke 
by 

b) Rudder deflection (zc) is defined as positive 
trailing edge left from the rear and is limited to 

Bee ce 15” 

c) Elevator deflection (n) is defined as positive 

trailing edge down (produces a nose down moment) and is 


iaimiited GO oS 
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d) Aileron deflection (&) is defined as positive 


With trailing edge of the right aileron un (produces a 
positive (clockwise from the mea) roll) and 1s J9m7sed to 
(EN 


Bier< 6 < 15° 


6) 


Se bette normzoemras stabilizer deflection (ny) ie 
positive trailing edge down. 
if vs 
f) Right horizontal stabilizer deflection (no) is 
Meseeive trailing edge down. 


fee ie 6 


B. ABRODYNAMIC COcFFICISNT BUILD UP 
Aereemnmamic forces and moments are dsuilt up from coeff 


Sm@emencata if Suvandard fashion using the following relations: 


a> GC. aps te Ug ac = Spee S 
aS fy ; = 
My = Oxy ase uy = GC, 450 vy = C, ol esis 


i) Beeanition of Coes 
Hach coefficient is partitioned into two parts; a 
Sueur c verm and a dynamic term as follows: 


Ca = 6 ne 


DD Ysp nv 
a CG. ate 
- Lon LoyN 
a eG 
LOY sp Lyn 
= Co + O 
ST DYN 
m= ¢ + ¢ 
1 Lep Lyn 


ao 








_ "ot "DYN 


Be GCoOocrimecrent Buald Up 
The static terms are constructed from coefficient 
Gata stored in tables (Appendix A) as functions of two or 


three variables (either <, 8, M, n, z,&}) as follows: 


é = f («,g,M) T . elevatar 
basic 

AG, oi = fh { <iehin) : 
ot 

Cy = f («,@,M) 
basic 

AC, (n) = eran) 
1, 

AC) (&) = f («,8,&) 
orl 

ACH ae ft Conner ) 
oak 

C =i (Gay B) 
basic 

ACy (zc) = f («,8,¢) 
Sl 

AC. (z) ae ( yen ) 
So 

C = f («,8,M) 

m ; 
basic 


AG (n) ae os exe) 
C = f («,B) 
AC (x) = f£ («,8,72) 


f cen) 


= 
2 
o— 

uN 
— 

I 
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C m=, 8) 


bas Cc 

Cc, (&) = f («,B,&) 
Suk 

Cy io) =f («,8,Z) 
5. 


teem stapmemecoctt@ecients are formed as: 


e 2-56 CAC ay) 
bem basic l 
C = +AC (n)+AC (&)+AcC (z) 
ee tate «= UST Den Dot 
Cc = VG +AC (z)+AC ( &) 
Yop = Ynasic)§=>0 Et top 
G =arC + AC (n) 
Mor "pasic "ST 
arn = Cf PeGnes) AG) (c) 
oT basic ST ST 
ci =a'6 + AC (E)+AC irae) 
Lom oo tnasic |= ott Lop 


Mme wien nhenmberms are Constructed from dynamic 
coefficient data stored in tables (Appendix A) as a function 
of angle of attack and Mach number f (2,M). They are 


Meee pte et My 
q q rp q 


Tee ue , GC 
ee L 


Cy fare ic Cy 


r p 


tee wtamiecwcOcthicirents are formed as: 


eG @ 
Ce ime fC; CQ) 
“DYN 2V, 7 aa 
Cc ° 
C oe. + Come Gg) 
'wyN  2V, °* De 
Cy a (Cj, ric, pt C,- 8) 
DYN 2V, *y p tg 





C = eee = + Cue -q) 


"DYN — 2Vp pee aa 
c = —— (ie Be Cc, rtc. p) 
"DYN 2V, OB r D 
Cy ge (Coe CP) 
DYN 2V,, r D 


The total coefficient build up is summarized below. 














C C +AC (n)+ s Cae °+C ) 
= NSD — = q 
L Lace Lon aV ao “q 
2 
aC +AC (n)+AC (&)+AC (Z )+ S (C~ «+0 q ) 
i Dee aT Dot Dog i eee 
. db a 
i ¢ mac. (c)+AC (= )+ (C, rtCy ptC.,. 8) 
t eee ST ton 2V Tp tp EB 
oe : 
C = C +AC (n)+ Cee q) 
5 ”" pasic om <Vin x q 
ts e. FAC (@RVtAC (Ej st (C, atc, r+C_ p) 
oe cacic ™ST ou Th a he a 
b 
OA = aC cemetery tAC, 8 (e)+ Cee sae oy, 
: basic ist i Ist <Vin +r 0 


GC. MISSILE LINEARIZED BQUATIONS OF MOTION 
The linearized quations of motion used to design the 
missile autopilot are presented in this section. The equations 
Meee imearized about the cruise flight condition. 
1. Additional Assumptions for Linearized Equations 
The additional assumptions required for linearized 
equations are as follows: 
a) Derivatives are given in stability axes which 
eooeticeemtor Lhe missile low level cruise condition. 
b) Standard small perturbation assumptions apply. 


1) Perturbed angles are small. 


oe 





2) Products of perturbed variables are neglected. 
3) There is no coupling between longitudinal 
end lateral directional motions. 
c) Thrust derivatives are neglected. 
d) The steady state condition is chosen as low 
meet Straight and level cruise. 
2. Linearized Equations Summarized 
The linear equations of motion in state variable 
memmeore Shown in lables2-1 and 2-2. The following 


Seeriicients are defined as used in Table 2-2. 


eel 
= ZZ XZ 3 XX 


ee = I 
TeX XZ XZ 


[ 7 r 
, eR ee a 22 

eee 

Zid oN XZ x 





ee tonto OlmmLabi lity Derivatives 


The dimensional aerodynamic stability derivatives 
Meee an the equations of motion are defined by the following 
relations. 


a. Longitudinal Dimensional Stability Derivatives 


ferns 420, )/at 1/sec 
U So Diy Dag Ss 
X= Tag (Cy -Cy )/m ft/sec 
59 a 
i io eweog.cc. /m ft/sec” 
A fe n So 2 
Y 
Peeps o. +2C~ Va 1/sec 
U 59 Ly Lag 55 
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2 2 ot I d d 
(nae2q)o = (Ana4?ayof fa] tConae1to © (Nat T)9 0 (2 na421)9 
(“Td+ N)V (TE+°N) Vy ( Td+ N)¥  ( Tat N)YV 0 (12 Tae 
0 0 ? 0 | 0 0 
Ss, SS, . Ss, SS, SS, SS, 
: a SS Qa" I YK © ksoo8 ox 
(ZL) (soTqeTtseA peqingsed) suotyenby TRucTIOOITG [Tes19y ey] 
2-¢@ ATAVL 
0 0 0 | 0 0 
S59 27 SDpy 27 SD A 27 SD 27 SD p 
eee) 3 ele a a Ss b fl 
Z Wt W Autss=W- (24°? N) Wt W 7 Wt OW Z° Wt W 
i 7 
on 27-98 9 7 -P? P=? Pn P7-PP 9 
47 Po hutsa- ar mr, a, 
4 i 59 4008 - 0 x ly 


(LL) (seTaetatea peqangsed) suotqgenby ier Ceol 


PC) Cae 
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Be = “FggS(C, + Cy )/m ay ee 
SS 
Ze = Bey, A 55 ft/sec 
2, = =GgSceC, /2n gq ft/sec 
’ = g. ae m ft/sec 
<p, 7n~ “dss ye / 
= qQ_ vy if 
My = FecgSG(C, +2Cy )/T Use 1/ft sec 
U Ss 
feewoaesscc.. /I 1/sec* 
« 4gs Me ey 
_- ee 
Mz = AggSe Gn. /Al yy gs 1/sec 
_- —2 
a = GgqSe ng! *tyy" ss 1/sec 
(,5 M_ = cee /i 1/sec* 
Me on © IssP°'m ! Tyy 
Pee baveral Directional Dimensional Stability 
Derivatives 
ea a 
Te = Iss S6y ,/™ ft/sec 
- 7 £ 
Y Fgg5bC,, /2nU a4 ft/sec 
Dp 
a = Bee Cy (7M ss fy asec 
Y- =Ge-SC, /m fy ft/sec* 
& i icc Ye Mm tO th sec 
a 4 2 
t. a oso y,’" 1, ity sec 
on, 2 
La = IggSPCy 0! Tx 1/sec 
oo Z 
L, = GagSb 1! *t xx’ ss i/see 
ns) sp2 eae 1/sec 
ie 1gs 1... oO 


> 








— he ae IggSPCq Tex 
Jr Le = Fg8b0q, /Ty 

N, > IggSPCn Tz 

Se :% 
i oss". ge, sé 
"Dp 

cc FggSb™ ¢ a /*t az" SS 
Ne oe ‘i IggSPOn Taz 
Mir XN = IgsShC h. 2 


Trim tHgquations 


eee 
Vecce 
eee 
1/sec 
1/sec 
eee 


1/sec 


Pomecurorehueand Level flight at constant velocity 


the following variables are all zero: 


=Z=0 QggrF 


Weg Sog"Sag 


The 


Tsint 4 


Tcos« 


coum 


ggg gr? 


95 


=0 


trim equations reduce to: 


nor 
O exe 


= mg-(C +C_ 


“3S 


ae 


| 
“-——™ 
CQ 


Dee Dees 


p Ngg)4g5 
O a YN 


7 *ss*On "ss)4ss5° 
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III. DEVELOPMENT AND SIMULATION OF THE MISSILE AUTOPILOT 


The missile is assumed to be equipped with an autopilot 
capable of providing closed loop control which consists of 
oe, Ollowing: 

a) normal acceleration (n) 

b) lateral acceleration (n) 

c) bank angle (6). 
ieeerdadition, an outer closed loop control, which serves 
as the outer autopilot, is as follows: 

meieeeetcitude (H) 

b) vertical flight path angle (y) 

These two outer loops are employed as required during certain 
phases of the attack mission. 

It is assumed that the missile possesses accelerometer 
sensors in the y and z body axes located at the missile 
eee GeOn gravity, rate gyros and rate integrating gyros 
about the x, y, and z axes, and a radar altimeter. Sensor 


noise is neglected in this simulation. 


A. AUTOPILOT INNER LOOP REQUIREMENTS AND DESIGNS 


1. Description of the Autopilot Inner Loops 


The autopilot is assumed to employ three inner loops 


as depicted in Table 3-1. 


ou 





LOOP 


1 Normal 
Acceleration 
Command 


2 Bank 


Angle 
Command 


Be burn 


Coordinator 


TABLE 3-1 


Autopilot Inner Loop Description 


COMMANDED CONTROL CONSTRAINTS 
VARIABLE 
n, (normal n 45° Se eae 
acceleration) Elevator 2a Ml Sa: 
e E 15> £ €£$ 15° 
(Bank Angle) Aileron > S$ 60° 


p S$ 50,100,200°/sec 


na (Lateral c St a 


Acceleration) Rudder ‘ae 


eeem@etes in Certain mission phases only 
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2. Normal Acceleration Command System 
The purpose of normal acceleration command (NAC) 
system is to provide a vertical load factor Va) response to 
a commanded load factor ie An accelerometer is used to 
provide the primary feedback anda rate gyro is used to 
provide inner loop pitch rate feedback for improved damping. 
a. Block Diagram 

The block diagram of the normal acceleration 
command autopilot is shown in Figure 3-1. Limiter a allows 
the commanded acceleration to a range from -2.0 to +4.0 g's. 
Limiter b allows the elevator control to vary from -15 
degrees to +15 degrees. 

b. Assumptions 

1) The rate gyro and accelerometer dynainic lags 
are negligible. 

2) The accelerometer is mounted at the center 
of gravity. Therefore the moment arm is zero and the 
feedback load factor is totally one ce=0)) 

3) The fix servo can be represented by a first 
order lag. 

4) The missile dynamics can be represented by 
a short period approximation. 

5) Commanded normal acceleration is limited from 
awe bo +4.0 o's. 

cay Design 
The following transfer function was developed 


from equation (11): 
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a(s) _ VaMntenMe)s + (Ma%y-ZaMy) 
n{s) ~ y Se cee.) cet GzeM -M.) 
T qx © wed. 
Vo Ver 


From this transfer function and equation (10a) 


V 
n= ee (Gas) 
& 


the following § to n, transfer function was obtained. 





2 
i = + (3 West (MZ 7} 
N(8) Vp -2 48° + (HZ tZ ple 8+ (Ma Zp -Zetty 
Q(s VaM +Z Me )st(M Z -Z MS) 
(s) g (Vy non & am ee a 


Peis tomemaitaemheriguics listed in Appendix A 
forethisg cruise missile at 0.75 Mach number, the anvove transfer 


functions become 





Q(s) _ -26.99(st1.009) 2,4 Nz'8) _ -0.08(s*+0.1985s-328.9) 
n(s) s°+1.2898+27.13 Q(s) (s+1.009) 


WIM ee ChaAraecyeristic polynomial of 
the Q(s)/n(s) transfer function yields the natural frequency — 
of 0.829 Hz and a damping ratio of 0.1237. The resulting 
block diagram for the normal acceleration command inner 
loop autopilot for analysis purposes is shown in Figure 3-2. 
(1) Root Locus Evaluation. 


(a) Pitch Rate Loop. The transfer function 


for the pitch rate loop is: 
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Ae 





Cie COM 
(s) (s?t41.29s*+(78.7+1080KR1)st+(1084+1080KR1) 





eo 


Appendix C contains the root locus plot of 
the pitch rate loop for ranging KR1. A KR1 of 0.28 is 
eieeen which yields the following characteristics: 


eae 5 OUU 23. SL 4. 5 


» = 55.53° 
r Oe 
ome |.06 
n a 


(b) Normal Acceleration Command Loop. 
ieeeecea hh! of 0.26 in the pitch rate loop and locating the 
acceleration at the missile center of mass (C=0), the normal 


acceleration command loop becomes: 





Be ee oes. 320) 
s(s°+41.298~+381.18-1386.4) 





Za =- 





N(s) Be GA (sO 23-329) 





ne (s) pGo 921. 2954381 . 15-1386. 4)-86. 4KA1 (8°+0. 28-329) 
C 


Hopendix © Contains the root locus plot of 
the normal acceleration command system for varying KA1l. A 
meior O05 as chosen which yields the following charac- 
teristics: 


eee eos 23.317 j 


b= 30.85° 
me 
el = 0.615H, 


423 








N (gs) 5411.29 





eauivecmeni Of O28 and a KAI of 0.05 the 


system transfer function relating os to a aes 
c 


N, (s) Wee 0k 25-329) 


34376.78 8©-1387.268-1421.28 


The normal acceleration command autopilot 
MeaneezerO SucCady Stave error for a step input as itis a 
type 1 system. The system response to a 1g commanded load 
factor is shown in Figure 3-3. 
3. Bank Angle Command System 
The purpose of the bank angle command (BAC) auto- 
pilot is to command missile bank angle. Limits are applied 
to commanded roll rate and in certain mission phases 
commanded bank angle. A rate integrating gyro feeds back 
bank angle and a roll rate gyro provides feedback as a 
roll rate damper. 
a. Block Diagram 
the bank anele command autopilot is presented 
PorOcmmatagram in igure 3-4. Limiter a allows the 
commanded bank angle to vary from 0 to 60 degrees and is 


aetive in certain flight phases. Limiter b limits the roll 


rate to either 50, 100, or 200 degrees/second. Limiter 


ec limits aileron deflections to range from -15 to +15 
degrees. 
b. Assumptions 
1) The rate gyro and displacement dynamic lags 


are negligible. 
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2) The fin servo can be represented by a first 
Omacr lag. 
3) The missile dynamics can be represented by 
a roll mode approximation. 
4) There are limits on commanded bank angle and 
commanded roll rate as noted above. 
e. Design 


Hquation (12) can be written to find the 


B(s) o(s) and W(s) 


meperal transfer functions ace rs)” aC 


3 


where 6 =Sor ¢ By allowing a roll approximation the 
emoter fUNCLiON relating aileron and bank angle becomes: 
>(s) Lo 
s s(s-L 
é (s) (s-L.,) 








iemrOoL ay the Origin indicates the missile does not have 
Bank angle stability without autopilot stabilization. 


mereuation Of the transfer function coefficients results 


im the transfer functions o(s) . euro and p(s) = 
ets) ast 6. 2) S 

oS The maximum missile roll rate can be approximated 
st6.2 

be nax (116) (15) 
ee ~P = Sana = +" = 280.6 degrees/second. 

max , 
L. one 


The block diagram of the bank angle command 


system is: 


G e KD 7 | KR2 So Eon I at g 
= = s+Z0 st6.2 : 


Roll Rate Loop 








Gomi ceremocic rvaluation. 
Ceeroemenate Loop. The transfer function 


for the roll rate Loop is: 


P(s) _ 4640 KR2 
as) Geet 34 6. 0KR2 





Appendix G contains the root locus plot 
Sueeme roll rave loop for varying KR2. A gain KR2 of 0.10 


yields the following characteristics. 


eee Ot 132.36 ~6s 
b = 59.96° 
aero O 
w= 46<47H, 


(b) Bank Angle Command Loop. With a KR2 


meee! 0 une bank angle command loop becomes: 






LO4KD 
6262712) 






mene bramanrer function is: 


ds) _ 464KD 
dCs) Be 72st LOLKD 
The root locus plot is shown in Appendix ¢ 


for varying KD. Choosing a KD gain of 10.8 yields the 


meollowine characteristics: 
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s = -9.98 + 9.58; 
o = 45.056° 
e712 

uy = 2.19 Hy 


Using gains of 0.10 for the inner loop and 
10.8 for the outer loop of the bank angle command system, 


the system transfer function G(s) becomes 


mc = Onl 2 
o (s) oe Oe 71252 5011.2 


ana the system is a type 1 system for a step input. The 
transient reponse is demonstrated in Figures 3-5 and 3-6 
for a commanded bank angle of 60 degrees. 
eee Coordinator 
Mee turn coordinayor (TC) is a lateral autopilot 
which vrovides a body directional load factor, oa iwOlea 


eommanded load factor, n . Like the normal acceleration 
C 


command autopilot it employs an inner loop yaw rate damper. 
Micwemmere Loop uses Unity accelerometer feedback and has 
the Capability of moment arm feedback if the accelerometer 
is located at a different position than the missile center 
of mass. 

a. Block Diagram 

The block diagram representing the turn coordina- 

tor autopilot is shown in Figure 3-7. Limiter a limits 


the rudder deflection (cz) from -15 to +15 degrees. 
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Figure 3-6 BAC Autopilot Response to 60° Commanded Bank Angle (Altitude Response) 
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b. Assumptions 
1) The rate gyro and accelerometer dynamic 
lags are negligible. 
2) An accelerometer is mounted at the center 
of gravity. Therefore the moment arm is zero. G=0 
3) The fin servo can be represented by a first 
order lag. 
4) The missile dynamics can be represented by 
peed roll approximation. 
Cree Dosw on 
The duteh roll approximation involves side slip 
angle 8, and yaw angle ¥. Equation (12) with » neglected 


(assuming a small C, results in the following set of 





B 
Poland ONS : 
(s) 
I B\S 
Se ees =v) —— e 
ee c(s) |. C 
92 (s) 
~N | s -shi WAGs N 
a aay gc 
6 , . c(s) | | 
iio euOtEOwiiic Geanstem LmyclLiOns can be 
extracted: 
Q(s) : Y_st(YN,-VqN,-Y¥ N.) 
= a a 
C(s) Vins (YL, +V pl) st (Lo Nt V oN -Y Ng) 
Y(s) Vol st(Z Np Yen, ) 
— er 2 | — ~~ 
ace s{ Vins -(ZotVoN ) st (YAN +VaNg-¥ Ng) } 


De 





) Val s+ (YN g-¥QNg). 


Ris 
2 
C (s) Vins “Tg tVoN st (Yen tV oN g-¥ Ng) 





Using the equation (10b) in Chapter I the 


transfer function for yaw rate to n_ is 
Z 


N (s) Vp Ys +(YN-Y NL) st(¥Ng-YoN, ) 
R(s) g  VqN st+(¥.NQ-YQN_) 


Substituting coefficient data obtained from Appendix A 


results in the transfer functions as: 


R(s) _ 2-14.5(s+0.11) 


c(s) s“+0.2458t14.5 


N(s) Mom de 20 Oe 20.9) 


R (s) a0 <1 1 

The resulting turn coordinator block diagram 
for analysis purposes is given in Figure 3-8. 

CUMnG@re ecto ut yo lua tion. 
(a) Yaw Rate Loop. The loop transfer func- 
| tion is 

MGs > SOs 0 11) 

E(s) s?+40.258*+(24.3+580KR3) s+ (580+63.8KR3) 

Appendix CG shows the root locus plot of 

the yaw rate loop for varying KR3. Selecting a gain KR3 


Poon, @eesults in the following characteristics: 


Seen 2. 339) 
6 = 56.85° 
a Oe o./2 
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b 56.85 degrees 
(oemenecooed@nator Loop. With a KR3 of 


feo the turn coordinator loop becomes: 


2 
( ) eee es —0,01s=-20.9) 


N " 
y S aIn TG San REE RRR > 
C ee eee 956. 352605m52) 
ceembarcim@arereyTunCGLiOnN 2s: 
N (s) -81,.2KA2(s~-0.01s-20.9) 
N (s) aia = 0) 7S ee) 5 +0601. 248m 2KA2)s + 
ec 


1697KA2 
ihewrooty Locus plot for varying KA2 is shown 
in Appendix ©. Selecting a gain KA2 of 0.25 yields character- 


mrencs as follows: 


s = -2.206+ 2.8073 
> = 38.16° 
r = 0.618 

Ww = 0.5682Hz 


The turn coordinator system transfer function 


G(s) becomes: 


N(s) RP so 00122009) 
N (s) BP yQNeee 221 2Mee 160). 55467803 


C 
The system is type 1 and therefore has zero steady state 


error to a step input. Figure 2-5a shows the response 
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of the turn coordinator for a step aileron input. Note that 
an ny of zero results in a coordinated turn. 
e 


B. AUTOPILOT OUTER LOOP REQUIREMENTS AND DESIGNS 
1. Description 
The outer Loop auto pilot has two functional loops 
consisting of an altitude hold and a vertical flight path 
angle hold. 
2. Vertical Flight Path Angle System 
The purpose of the vertical flight path angle hold 
Meopeor famma command Loop is to provide missile reponse to 
a command of 8.6 degrees during a specified phase in the 


Mmemmmnal profile. 


a. Block Diagram 





b. Assumption 
1) Bank angle is zero. 


c. Design 


V 


Starting with equation (10a) ae “n(Q-«) eee leey, 
24 


g 
y(s _ 
the transfer function rte) = Tats) can be derived. Using 


wc titer Loop Gecten pmeyiously for load factor the vertical 


Plight path angle block diagram becomes: 


a 











-9.504(s~+0. 28-329) 


34+41.298°4371.68~+1 384843127 





wm |— 





ae transier function ——- 1S 
Cc 


Vis) _ -9.504K (s*+0.2s-329) 


Me) e+). 2947437. 687+ (1384-9. 504K )s*+(3127-1.9K |) s#3127K , 





(1) (Relot—heeus eee on’. Appendix C shows the 
meomeLocus plot Of the vertical flight path angle loop for 
varying gain a A gain Ky of 1.0 yirelds the following 
characteristics: 


Sieh oe 2.5 47 7 


> = 26.6° 
c = 0.4478 
De 0. 45H, 


The system transfer function G(s) with a 


Som Or 0 becomes: 


ys) _ -9.504(s*+0.2s-329) 
v, (3) BERN 7 oe 27). 53-4 3125. 1843127 


Figure 3-9 depicts the system response to 
a commanded vertical flight path angle of 8.6 degrees. 
Sumetitude Hold System 
The purpose of the altitude hold autopilot is to 
hold altitude at 50 feet during specified flight phases. 


a. Block Diagram 
K 


Ba teal VAS) H(s) H 
. c [Vm on Y (s) 
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be Assumptions 
1) Vertical rate of climb (or desent) 
subtends a small angle (yis less than 20°), therefore 
H = Vpsin y = Vinv. 
2) Bank angle equals zero. 
ec. Desicn 
Let altitude command and altitude be designated 
Ho and H respectfully. The design is based upon the 


following two relations: 


H = Ki(H ¢d) 
K, 
Yo = gt (ASH) 
T 


From the vertical flight path angle development, 


a Te 
N ot Vins 


Z 


sH(s) = Vas) (Divide both sides by a(s)) 


sH(s) _ V s) = Vy g 
Ns) I Ns) a Vins 


H(s) 


y(s) S 
Tien aumoudeo nold biock diagram Ddecomes: 


~9.504(s~-0.2s-329) 1| 4 


Bee aos eo 12s. 1st3127! |s 





The transfer function H(s)/H,(s) is: 
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Hg 29. 504KH(s“+0.2s- 329) 
7m 7 4 ) 2 


Tm 
s) $°+41.295°+371.68 941374. 587+ (3125.1-9. 504KH) 8 “+ 
nei ki)st3127Ki 
(1) Root Locus Evaluation. Appendix C contains 
eeecreooc focus plot of the altitude hold loop for a varying 
mene. Acain Ke of 0.3 results in characteristics: 


s = =-1.4045 + 2.4243 


b = 30.08° 
aa Ono cn 
as =e wea. 117, 


The altitude hold loop transfer function 


Mesh with KH of 0.3 is: 


H(s) _ 2.8512(s*+0.2s-329) 
7 5 4 2 


H,(s) mer. |.297s +371.6s 41374.5s mer 22. 256+ snl2omet 938.1 


and Pigure 3-10 shows the response to ay, of 50 feet. 
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IV. DEVELOPMENT AND SIMULATION OF MISSILE GUIDANCE SYSTEM 


ieee scewOr cals Study is to evaluate the effect of 
limited roll rate and various guidance laws on the terminal 
performance of a bank-to-turn cruise missile using propor- 
tional navigation in the terminal (attack) phase against a 
medium to large combatant ship. A baseline guidance and 
control scheme is designed and studied which uses a pop 
out maneuver, proportional navigation, and bank-to-turn 
maneuvering in the attack phase. The trajectory parameters 
and navigation constants are optimized for best performance 
against a moving ship utilizing blinking countermeasures. 
The missile is assumed to have limited roll performance 
Smemeoo bie use of differential stabilizer for roll control 
vice ailerons. The terminal performance of this baseline 
system is evaluated for varying roll rate limits, electronic 
eountermeasure blinking rates and burn through ranges. 

meocecoma control scheme 1s also designed and evaluated 
for the attack phase. This mode uses a sea skimming scheme 
in which a 50 foot altitude is maintained and proportional 
navigation is used in azimuth during ingress until a specified 
range is reached at which time proportional navigation is 


used in both azimuth and elevation to attack the target. 
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Meme MicoanLON DESCRIPTLON 
1. Baseline Guidance Scheme (Pov Out Maneuver) 
The mission is divided into four phases: 


a iggyseaucycus 


Oe 


© 


popup 


) 
) offset 
) 
) attack 


Qu 


The mission commences with target ship located 25,000 feet 
north of the missile and moving east at 21 knots. The 
Missile initial heading is north. 

ime teres piase Commences with the missile flying 
at an altitude of 50 feet on altitude hold and utilizing 
Proportional navigation in azimuth to home in on the 
target. Ata given range the offset phase is initiated with 
“mee borne right of approximately 12 degrees followed by 
Mumeonindnaeda bank angle of zero degrees. No proportional 
Perea vron 15 Used in this phase. At another given range 
&@ pop up maneuver is initiated which generates a low angle 
climb in preparation for the attack phase. During the pop up 
maneuver, proportional navigation is used to begin turning 
Mie Missile toward the target. Finally the missile enters 
Missdvuaee Nase ae a specaified’range where proportional 
Mev catien is Used mm both vertical and azgimuth planes. 
im this phase the bank angie is unlimited but roll rates are 
limited and used as one of the parameters in studying miss 


distances. The second parameter under study is the electronic 
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counter measure blinking rate frequency which is initiated 
at the start at the start of the attack phase and allowed to 
Pemain activated until a specified burn through range which 
is the third miss distance study parameter. Figure 4-1 is 
a geometric depiction of the missile flight path. 

2. Alternate Guidance Scheme (Sea Skimmer) 

The mission for this guidance scheme is divided into 
two phases: 

ae ingress 

be. attack 
This mission begins with the same initial scenario as the 
baseline scheme, however, the pop out maneuvers (turn and 
pop up phases) is eliminated. : 

The missile flies at 50 foot altitude during the 
ingress phase and utilizes proportional navigation in 
azimuth until a range of 3000 feet at which point the 
missile enters the attack phase and dives in on the target 
Using proportional navigation in both azimuth and elevation. 
This low altitude approach essentially resembles a sea 
Skimming effect. Figure 4-2 shows the geometric flight 
profile of this scheme. 

3. Electronic Counter Measures (CM) 

it is assumed that the target has the HOM capability 
of simultaneously shifting the aim point of the missile 75 
feet aft and 10 feet vertically. The simulation is tested 


Pec Glieorinkine frequencies of 0.5, 1.0, and 2.0 cycles 
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Figure 4-1 Baseline Guidance Scheme Flight Profile 
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Figure 4-2 Sea Skimmer Guidance Scheme Flight Profile 
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peeecoceconce Blinking is commenced during the attack phase 
and continues until jammer burn through occurs when the 
target signal power can be seen over jammer power. Beach 
Beme feet (as calculated using Hosington, {Ref. 3}) is 
considered a typical burn through range the missile is 
tested for burn through ranges from 3800 feet down to 
400 feet. 
4. Glint 

Mees Souey models glint with a gaussian function 
that gives the aim point a random shift effect with a 
Poeewemce Of 25> feet in the Y axis and 4 feet in both the 
Keana@ Z axes. Since glint is imposed as a fluctuation in 
target coordinates and target coordinates are used to 
Calculate the target line of sight angles then glint is 
modeled as inversely proporticnal to range as required. 
For ease in simulation glint is imposed only during the 
Pereen phase. Results for runs with glint included are 


presented in Chapter IV. 


B. SBEKER EQUATIONS AND SIMULATION 
A flow diagram of the seeker simulation is presented 
in Figure 4-3. 
io ine Cr coment Rates 
It is assumed that the missile seeker has the ability 


meneoense tarpet range and target line of sight rates in 


azimuth and elevation (6 and on, ) in the body coordinate 
Lap BL. 


system out to a range of 25,000 feet. 
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These line of sight rates are generated in simulation 
ieecwO atuep procedure. First, the line of sight rates 


) 


in earth referenced azimuth and elevation ere and o ~ 


co 


meomconstructed by using the following relations: 


¥ Lo 
a, = tan™! eee 
Tn Kn 
_ -1 -H 
een ee 
; ((Xp-Ky) “+ (Ya-YQ)*)® 
= : 
AZ *t : 
a7 = tno 
nr = U7 9 
oe = ° _9 
BL a 


Figures 4-4 and 4-5 depict the geometry represented by 
the above equations. 
meow vic Dodyelane of sight rates are calculated from 


the following transformation. 


2. = - On, Sin d+ Ong cosgcoso 

e a @ @ 9 a 

EL. 047.cosot+0,,cos sin® 

Onee the body rates have been simulated (6 ole ) 
AZy BL, 


bhey are transformed to earth horizontal and vertical rates 
through the transformation Ty defined by: 


oo. Syeoke -~¢ 
a; yh a a 


oO... Bae 
AZ ~ BL, hl cos? 
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Figure 4-4 Azimuth Line of Sight Description 





6 
T (shown negative) 


(y*4y~) 2 


Blevation Line of Sight Description 


Figure 4-5 
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This transformation includes the assumption thatQ@is small 
Baroughout the mission. The earth horizontal and vertical 


line of sight rates (o,, and Ong) are then used in the 


mepeopriaye fsuidance laws for proportional navigation. 
ee Vaiccile Dynamic Filters 
imocdoreconhilter Out missile dynamics from the 
seeker line of sight rates, a low pass filter, 10/(S+10), 


is used at the input to the seeker on 9 and wv. 


C. BASELINE GUIDANCE LAW DESIGN (POPOUT MANEUVER) 
A flow diagram of the baseline guidance scheme is 
saown in Figure 4-6. Switch positions are a function of 
mission phase. Hach phase of the mission is described 
below beginning with the attack phase. 
1. Attack Phase 
Proportional navigation is used in azimuth and 
elevation to construct guidance commands to the missile 
autopilot commencing at a range of 9100 feet in the manner 
described below. 
a. Proportional Navigation Constants and Guidance 
Commands 
ine carvueiominwental and yertical line of sight 
rates ies and a which are generated from the body rates 


and Gp, by the transformation T, are multiplied by 


B B 0 


a navigation constant and Val g to generate commanded lateral 


a, 
accelerations in the earth vertical and horizontal planes. 
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A navigation constant of 4.0 in each plane was found to 

yield best results. Values from 3.0 to 4.5 were evaluated 
with results of that evaluation shown in Chapter V. A one 
feomes is applied to the vertical commanded lateral accelera- 
tion to compensate for acceleration due to gravity. Thus 


two lateral acceleration commands (n, and ne ) are computed 


e € 
in the earth horizontal and vertical planes. 
n = xX a 5 
te AZ g AZ Ge) 
Vin " 
‘a tte, 
C g 


ck 


These commanded lateral accelerations are subjected to a 
further transformation (Ty) in order to generate commanded 
bank angie and normal load factor which are required by the 


Meese autopilot. 





-1 ee 
o a Tea 
n 
Z, 
Cc 
n, = ne, cos otny, Sing 
e e e 


Peeure 2>/ depicts the geometric relationship described by 
the above two equations. Commanded bank angle (>,) and 


normal load factor (n ) are used as inputs to the appropriate 
C 


autopilot outer command loops to provide closed loop bank 


angle and load factor responses. 
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b. Limits on Guidance Commands 
In the attack phase, the following limits are 
applied to commanded variables: 


a2 n Ss 


Missile performance is evaluated for maximum roll rates of 
memeero0 and 200 degrees per second. 

Cm banksanele Anomoly 

InStne atuack mode there is no limit on 

commanded bank angle. Thus a computational scheme is 
memercd LO insure that the missile always rolls the 
shortest way to a new commanded bank angle. This is accon- 
@eoted by continuously calculating the quantity: 


a= iG at) = ye 


If A(t) is less than 180 degrees then > (+) i Sera tL Oma 

bank angle autopilot. If A(t) is greater than 180 degrees 

then a modified bank angle command (>, eo schimnor ine. 4 ulto 
M 


Pulot in accordance with the following logic: 


d a ,+360° for > <0 
CQ 


b = 6 -360° for 6 >0 
¢ c 


d. Switch Positions 
The switches shown in figure 4-6 are in the 
mou LOwWinNg Positions for the attack phase: 


SW father position 


fy 














| 
{ 
{ 


SWe Peostulon < 
SW 3 Poorer on | 
SW4 Closed 
e. Attack Phase Summary 
In the attack phase proportional navigation is 
used in both azimuth and elevation to generate bank angle 
and normal load factor commands to the missile autopilot. 
The phase is commenced at a range of 9100 feet with missile 
Ofiset in heading from the target at an altitude of approxi- 
mately 250 feet. Roll rate is limited to ORE Bank angle 
monmou limited. Normal load factor is limited from -2 to 
eee S 5 
oeeeress Phase 
Miewiaieress onace commences with the beginning of 
femerolLem Simulation when the missile is 25,000 feet from 
the target and initially heading north. The missile's Mach 
number is 0.75 and the ship's velocity is 21 knots to the 
east. 
a. Altitude Hold Methodology and Guidance Commands 
The missile is required to maintain altitude of 
50 feet above sea level_.during this phase. It accomplishes 
this by — ee present altitude with 50 feet and multiplying 


the error by —~ to establish a flight path angle command. 
T 


om 
angle command is compared with the present flight path angle 


is assigned a value of 0.3. This calculated flight path 
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V 
wiGmoie Gimterciece Om error is multiplied by K ,and = 


to give an elevation line of sight rate, 0 or - A value of 


emis assigned for Kye ats ne MoncOlpuLGed asumaollows: 


C 
Ven 
his = Key a oe) | 


The lateral acceleration command Ny ioe e Wino a 
c 


exactly as described in the attack phase with bank angle and 
meme rates limited. 
The second transformation for the ingress phase 


(T..) remains the same form as the attack phase. 
N p 


4 


b = tan on 
C d 
ene 
n 
Z, 
Cc 
n = Ny cos én, Sind 
“@ eC C 


b. Limits on Guidance Commands 
PiPiriememevess phase the  rollowing limits are 


applied to commanded variables: 


< < 
-2 =e ae 
C 
lo | <lo! 
joo | <lhax | 


The maximum bank angle allowed is 60 degrees and roll rate 
as fimited to 50,100 or 200 degrees per second. 
C. owWitch Positions 
The switch positions as shown in Figure 4-6 are 
as follows during the ingress phase: 
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SW1, SW2, SW3 Position 1 
SW Closed 
3. Offset Phase 
The missile initiates the offset phase at 17,800 feet. 
The offset phase is subdivided into two sub-phases -- turn 
and level. 
a. Turn and Level Sub-Phases 
The missile maintains an altitude of 50 feet 
jeouehout the offset phase in the same manner as described 
during the ingress phase. Thus 
a ees Ge) es 
g 
The missile performs the coordinated turn by 


using a commanded bank angle and computing a normal load 


mee wor as follows: 





Pe waae 
C - 
z 
Ce 
al = 
mc Cos! ip 


At a range of 15800 feet from the target a turn of approxi- 
ieee \zncerrcoos is accomplished. The missile then performs 
a straight and level maneuver by commanding zero bank angle 


along with the same normal load factor as follows: 


a FAG 
D eae 0 


SQ 








b. Limits on Guidance Commands 
Poe ewornse up ohase normal load factor is limited 


aoe OlLlLows: 


Commanded bank angle is fixed at 60 degrees during the turn 
sub-phase and 0 degrees during the level sub-phase. 
Ce Brronwrost t1LOns 
The switch position, as shown in Figure 4-6, are 
fomeo!lOws during the offset phase: 


1) Turn sub-phase 


SW1, SW2 Posmuron .| 
SW 3 hocastaom 3 
SW4 Open 


2) Level sub-phase 


eee Rosier On le 
SW 3 Posi tmen 
SW4 Open 


wooo Up Phase 


The missile enters the pop up phase when the range 
becomes 10,200 feet. The purpose of this phase is to 
Popup udemmmsscile to an altitude of approximately 250 feet 
so that a dive into the target will give acceptable deck 
penetration. 

ae cocebedamritoht Path Angle Methodology and 

Guidance Commands 
The Missile achieves the pop up effect by 


commanding a vertical flight path angle of 8.6 degrees. 
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This commanded climb angle is compared with the actual y and 
mie erroreis multiplied by KY and 7 to give an elevation 
Mme Of cicht rate. A value of 1.0 is assigned for Ky 


iechis simulation. Thus ne is computed as follows: 
e 


T 
eee ‘ zs 
ZL Y ( climb 

Cc g 


Y) + 14 


The Lateral acceleration commanded Mh is compuced 
é 
exactly as described in the ingress phase. 


The second transformation Ty iS eclS Ode Ono. : 


a _ 
O = ayer Yo 
e -— 
e 
n, = ne cos otny sind 
© e c 


Vi ies io cc uUOnSs tart rolling into the 
target using proportional navigation during the pop up 
phase and thus aids in the attack phase which follows the 
Meme up phase. 

b. Limits on Guidance Commands 
Polen eOpeUpmpicse the following Limits are 


applied to commanded variables: 


<a < 
a 23 Zs 
a 
Poo a aa | 
lo | & Jo | 
. max 
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The maximum bank angle is limited to 60 degrees 
and the maximum roll rate is limited to 50, 100 or 200 degrees 
per second. . 

eo paviveh Positions 
The switch positions as shown in Figure 4-6 are 


as follows during the pop up phase: 


SW1 Pocwmrion 2 
SW2, SW3 Eooiemen | 
SW4 Closed 


5. Guidance Summary Baseline) 


ioe tOllowine Watrix indicates the horizontal and 


fmeeuieal Controls during each phase of missile flight. 


PHASE AZIMUTH CONTROL ELEVATION CONTROL 
ingress PN* abate Oe ou cd 
Offset ®o Command Altitude Hold 
Popup PN Flight Path Angle Hold 
Attack PN PN 


*“Proportional Navigation 


D. ALTERNATE GUIDANCE LAW DESIGN SBA-SKIMMER 
The sea skimmer scheme differs only from the baseline 
guidance scheme in that the offset and popup phases are 
Mma te@s Im csosence this scheme is composed of an ingress 
phase and an attack phase. 
I Gaeress Phase 
The ingress phase commences with problem initialization 


which locates the missile 25,000 feet from the target with 
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| 
| 
: 





mlenmbical initadl Conditions as the baseline guidance scheme. 


The vertical load factor is computed as 


Ver 
Ny = Key (Y = Y) + 1 
C s 





C 


The lateral acceleration command is computed as 





The second transformation (Ty) io eas ep Olle ss 


d = tan”! “c 
C n 
Z 
C 
and n, = Ry cos On, sing 
C C C 


Limits on guidance commands and switch positions 
Meio ude rcss phase are identical to the ingress phase 
meee baseline scheme. 

e- Attack Phase 

The attack phase commences at 3000 feet range and 
uses proportional navigation in both elevation and azimuth. 
The missile dives into the target from 50 foot altitude. 

The bank angle command and normal load factor are 
computed exactly as in the baseline scheme attack phase 


eae are: 


SA 





=——eeosorn.. Sin 
n, 7 p y Si Wy 
Cc Cc Cc 


[fitusmeonecuraance commands and switch positions 
are also identical to the attack phase in the baseline 


guidance scheme. 


3. Guidance Summary (Sea Skimmer) 


iMemioluowine Watrax shows the horizontal and vertical 


emeetols Guring cach phase of the sea skimmer guidance 


aomeme . 

PHASE AZIMUTH CONTROL VERTICAL CONTROL 
Hea co THESIS PN HAL Ier. ives ali 6. 
Attack PN PN 


So 





V. PROBLEM SIMULATION, RESULTS, CONCLUSIONS AND RECOMMENDATIONS 


This chapter describes the computer simulation special 
language used and presents the results of test runs in which 
mom rate, BOM blinking rate and burn through range are all 
varied in both the baseline guidance scheme and the sea 
skimming guidance scheme. From these results optimum para- 
meters are chosen to yield minimum miss distance for each 
guidance scheme. Miss distance ranges are also presented 


for the baseline scheme with Glint imposed. 


A. CSMP SIMULATION 

The Continuous System Modeling Program III (CSMP III), 
an IBM developed system, is the modeling technique chosen 
Pees suudy. OsMP, as discussed in Speclchart and 
has six selectable integration methods (fourth order Runge- 
Kutta is used for this simulation), and has multiple special 
Mmuerton codes that can solve a large variety of control 
engineering problems. 

The special functions used for this simulation are: 

a. FUNGEN - This function extracts the value of a 
function of one variable from a stored table in memory. It 
is used for aerodynamic coefficients. 

b. TWOVAR - This function extracts the value of a 
Function of two variables from a stored table in memory. It 


is used for aerodynamic coefficients. 


86 





e. INSW - Y = INSW (X1, X2, X3) assigns Y = X2 for 

x, <0, Y = 3 for X, 2 0. It is used to accomplish ECM 
Geraking action. 

Q@. GAUSS - This function is a random number generator 
Mmeenormial distribution. it is used to realize glint. 

a werbeb = this funetion simulates a first order 
ieee; (tS*+1) Filter. It is used in autopilot and seeker 
SBeweacions as low pass filters. 

f. DEBUG - This is a computer program aid that evaluates 
all output parameters and variables at requested times. 

mieaeeartion vo the special functions, trigonometric, 
integral, derivative and square root functions are used to 
perform necessary steps with one invocation. CSMP also has 
the ability to perform calculations in either a sort section, 
where the order of computation is not considered, or a 


Motte Seemwron, where ordered sequential fortran statements 





Seemreguarcd, Ihe program is structured as shows below: 
MISSIL& SIMULATION PROGRAM 
INITIAL 


CONSTANTS 
APRODYNAMIC COEFFICILNT TABLES 


DYNAMIC 
NOSORT 


MISSION PHASt LOGIC 

GUIDANCE BQUATIONS 

ECM (BLINKING) EQUATIONS 
SORT 


6 DOF MISSILE DYNAMIC EQUATIONS 


Oo 








INNER LOOP AUTOPILOTS 
SELKER EQUATIONS 


TERMINAL 

OUTPUT GRAPHS 

The CSMP initially sets up the two and three variable 
Seewticieénmus in tabular form in memory. The run commences 
With the missile in steady state flight with steady state 
variables established. These steady state variables are 
used to enter the tables in memory and extract coefficients 
which are used to calculate aerodynamic derivatives, forces, 
fememus, and Ouneér relations. From the steady state flight 
condition small perturbations are made and new table entries 


meoemeed. Ihis process is repeated for all phases of flight. 


B. RESULTS 

Both the baseline guidance scheme and the sea skim 
guidance scheme were tested under the same conditions and 
quantitative results obtained. All scenarios started from 
eo OlLVowrie initial conditions: 

(Peed eet position - 25,000 feet north of the missile. 

2. Target speed - east at 21 knots. 

3. Missile speed - north at 0.75 Mach number 

4. Missile altitude - 50 feet above sea level. 
The simulation was run for 31 seconds. Miss distance data 
was obtained for varying roll rate limits, ECM blinking 
Paecwhemcmes, ama burn through ranges. The following matrix 


of parameters was tested: 
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Burn through range (feet) 
2000 
deo 0 
1 2100 
800 
400 
Roll rate limit (degrees per second) 
200 
ee 
50 
ECM blinking frequency (Hertz) 
Z 
| 
5 
Piewbesulting Miss distances are shown in fables 5-1 
and 5-2 for all runs. Figure 5-1 through figure 5-6 are 
the graphical results of simulation runs. Each run is 
presented in 6 consecutive figures a through f and depict 
the following. 
a. Commanded roll rate (PCMD) and roll rate (P) vs 
time. 
b. Commanded bank angle (PHICMD) and bank angle (PHT) 
vs time. 


c. Range, X Range, Y Range and Z Range vs time. 


Peeteiweand righ, stabilizers (LSTAB and RSTAB) vs 


e. Rudder position vs time. 
f. Map of YMISL& and HMISLE vs XMISLE. 
The following matrix indicates the test conditions used 


to generate each figure: 
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ECM 
mer nking 
Rate 


re, cps 


im eps 


O cps 


TABLE 5-1 


Baseline Guidance and Control Scheme 


ont 
Rate 
arn 


DO 
100 


200 


50 
a0 
200 


De 
100 


200 


50 
100 


200 


Miss Distance (feet) 


400 


Voy 
HP 
Zi 


24 
12.6 
D4 


2D 
ed 
13.3 


BS, 
44 
ae 


Burnthrough Ranges 


800 


Sef - 


Tio 


30 


ied 


20 


eo 


79 


90 


1200 


16 


a7 


sietlers 
ae 


a) 


1500 


fi 
oAo> 
ad 


iyo 


tal 


ae 


eno 


Aer 


2000 


O25 
ie emne 
Oem 








CM 
Blinking 
Rates 


Sea-Skimmer Guidance and Control Scheme 


roll 
Rate 
ama, & 


50 
100 
200 


50 
100 


200 


50 
100 
200 


50 
100 


200 


400 


43 
We 


44.64 


Be 
De 
SD 


34 
Be. 
43. 


TABLE 5-2 


Miss Distance (feet) 


O- 


nN 


800 


ers 


61 


aie 


D4 
De 
63 


oO 


Ae 
ee 


Burnthrough Ranges 


a 2210, 


bo 


91 


1200 


37. 


61 


Dre 


38. 


20 
26 


61 


ae 


62 


1500 


Dic 


61 


ie? 


Coe 


AO 


26 


Suile 


61 


28) 


On 


2000 


toe 


iS 


B09 


pe 











BAS& GUIDANCE SCHEME 
5-1 (a-f) NO ECM, 100°/see roll rate 
5-2 (a-f) ECM (1Hz), 100°/sec roll rate 
5-3 (a-f) ECM (1Hz), 50°/sec roll rate 
5-4 (a-f) ECM (1Hz), 200°/sec roll rate 
SHA SKIMMER SCHEME 
5-5 (a-f) NO ECM, 100°/sec roll rate 
peo (a-f) ECM (1Hz), 100°/sec roll rate 
[eee malysic Of Bascline Scheme Results 
a. Varying Burn Through Range 
Miss distances generally decrease as the burn 
GarOUgh range increases. For runs beyond 1500 feet burn 
through ranges, the miss distances are close enough to be 
considered as hits since most of the miss distance is in the 
Y component of range which has the largest offset (75 feet) 
applied during s0CM. 
b. Varying ECM Blinking Frequency 
ine Wiss distamece decreases as blinkiing rate 
increases for the same burn through range. Miss distances 
for burn through ranges beyond 1200 feet are varticularly 
peeGgs ror dll three roll rates. 
c. Varying Roll Rate Limit 
Varying roll rate limit results in the largest 


Variation in miss distance. For 200 degree per second 


Gora comma, and 0.5 Hertz BCM blinking frequency the 
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Sea Skimmer Response for ECM 1 CPS, Roll Rate Limit 100 deg/sec 


Figure 5-6c 
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(termediate s Temains high for ail burn through ranges out 
to 2000 feet. But, for the remaining blinking frequencies 
peo DUG Lorough ranges the miss distances improve signifi- 
Satuly fOr Inercasing roll rate limits. A miss distance of 
7.9 feet for a burn through range of 800 feet and 2 Hertz 
Pmmmcams frequency is excelient. 

2. Analysis of Sea Skimmer Scheme Results 

The results for the sea skimmer scheme are quite 
Gitferent than the baseline scheme. 

Burn through range seems to be the single most 
Teron vaedabives wMiss distances for a burn through 
range of 1500 feet are marginally satisfactory for all roll 
rave 2imits and HOM blinking rates. However, ata burn 
feedmranse Of <000 feet, miss distances for varying roll 
Peeomeatmbts and GCM blinking frequencies are extremely low 
and considered quite good. 

Varying the ECM blinking rate produces little effect 
on the miss distances. However, a unique roll rate limit 
of 100 degrees per second produces miss distances double the 
miss distances for roll rate limits of 50 or 200 degrees 
mer second. 

3. Analysis of Glint Simulation Results 

Figure 5-7 shows the random introduction of glint 
Bimodewtameet POSLtioOn and line of sight rate. The miss 
Gistance for the run in this figure is 15 feet, but analysis 
from other runs show that glint increases the miss distance 


poe co 50 "percent. 
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4. Proportional Navigation Ratio Selection 


The following matrix shows the result of miss distance 


mieeeno BOM environMent for varying proportional navigation 


constant (PN). 


PN Miss Distance (feet) 
3 0 © 
345 6.9 
4 4e4 
a) Dia 


Although a proportional navigation constant (PN) 
of 4.5 results in the smallest miss distance excessive over 
MiesueLoeobeD inpuLs results in this proportional navigation 
fwmebant as being an unsatisfactory choice. A PN of 4.0 was 
chosen based upon satisfactory system responses and the 


small miss distance achieved. 


C. CONCLUSIONS AND RECOMMENDATIONS 

Based upon the results of this study it can be concluded 
that miss distance is a function of roll rate limit, ECM 
blinking rate, burn through range, and the guidance scheme 
ehesen. the guidance scheme chosen is most significant 
Of these factors. 

For the baseline guidance and control scheme, roll rate 
PiiieLs ome Most Significant variable affecting miss 


distance. One exception to this is the 0.5 Hertz blinking 


Ie 
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Berselime Scheme - Koll wate limit = <00 degrees per 
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wie LOllomie adcitional studies are recommendec: 
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Ds @fiation of the cirection and magnituce of the 
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APPENDIX B 
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BASELINE GUIDANCE LAW SIMULATION 
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